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ABSTRACT 

I present details of the variations of several hundred red giant stars on time 
scales of a few hours to a few days from Hubble Space Telescope (HST) observa- 
tions of a low-extinction galactic bulge sample from an intensive seven day cam- 
paign. Variations in the red giants are shown to be a strong function of position 
within the color-magnitude diagram (CMD) in accord with general expectations 
from theory. Amplitudes are greater for stars with larger radii, whether this 
results from higher luminosity at the same effective temperature or lower tem- 
perature at a fixed apparent magnitude. Likewise, characteristic time scales for 
the variations increase to the upper right in a CMD as does the ratio of ampli- 
tudes measured at 606 nm compared to 814 nm. Characteristic variation time 
scales are well matched by low-order radial pulsation modes. The effective sam- 
ple discussed here extends from about two magnitudes above the bulge turnoff at 
which red giant radii are ~7 R/R Q at 5,000 °K with typical amplitudes of ~0.5 
mmag to ~40 R/R & at 4,000 °K with amplitudes of ~3.5 mmag. Variability 
characteristics are quite similar at any given position in the CMD, and at lev- 
els in the CMD where oscillations are easily detected nearly all red giants show 
such. If these variations represent oscillations with sufficient lifetimes to derive 
accurate mode frequencies more extensive observations, e.g. as should soon be 
provided by the Kepler Mission, would provide a rich asteroseismic return. 

Subject headings: stars: late-type — stars: oscillations — techniques: photomet- 
ric 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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Introduction 



Red giants are attractive targets for the study of stellar oscillations. From the observa- 
tional perspective, the am plitude of stochastically-dri ven p-mode oscillations is expected to 
rise with increasing L/M ( IKjeldsen &: Bedding! Il995l ). Assuming that coherent oscillations 
exist with sufficiently long lifetimes to allow studies analogous to those from solar oscillations, 
asteroseismology of red giants promises probes of great interest for testing the theory of these 
highly evolved stars with complicated interior structures. For the red giant phase, where 
evolutionary tracks from different progenitor masses converge in a CMD, asteroseismology 
can, in principle, provide valuable const raints on the intrinsic parameters of individual stars 
rtstello et alJl200a iKallineer etaDhoOsh . 



That red giants show characteristic variations is well known. Less clear (and the subject 
of increasing study) is what the physical nature of these variations are. A central question 
for variations in any group of stars is whether these follow from normal modes of oscillation 
either driven by nonlinear feedbacks such as in RR Lyrae stars, or from stochastic driving as 
in the Sun. The common alte rnative to oscil l ation s would b e variations arising from g iant 
cell convection as suggested by iSchwarzschildl (119751 ) , see also iDziembowski et al.l (120011 ) and 
Ludwig! (120061 ). in which case their utility as probes of interior structure likely do not exist 
as with true oscillations. A secondary question is whether the variations, if oscillations, have 
lifetimes long enough to allow the accurate determination of frequencies required to support 
asteroseismic interpretations. This paper, while not resolving these questions, will more 
firmly establish the overall characteristics of red giant variations by utilizing a serendipitous 
data set that provides detections on several hundred stars within a narrow confine of the 
CMD. 



After a decade in which the promise was recognized (IBrown fc Gillilandl Il994l ). but 
results did not materialize despite concerted efforts, asteroseismology has provided solid and 
exciting results in rece nt years manifested by preci se measurements of solar-like oscillations in 
a number of stars. See lBedding fc Kjeldsenl (120071 ) for a recent general review and references. 
For red giants in particular a number of radial velocity studies requiring dedicated use of 
at least moderately large tele scopes and excellen t spect rographs have provided evidence 



for oscillations in red giants. Hatzes &: Cochran ( 1994al b) f ound 14-54 m s 1 variations 



in a Boo, following up on the ISmith. McMillan &: Merlind (119871 ) detection, and (3 Oph 
respectively using observations over typically eight consecutive nights on the 2.1m McDonald 
telescope. These velocity amplitudes would c orrespond to photometric variations of 0.5-2.0 
mmag using the IKjeldsen fc Bedding! (119951 ) scaling relation. With observ ations oyer 12 
nights and precisions better than those successful for the Kl III star a Boo, iHornerl (119961 ) 
failed to detect evidence of oscillations on any of four G8 III to K2 III stars studied. More 
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recently lFrandsen et al.l (120021 ) detected solar-like oscillations in the G7 III star £ Hya peaking 
at about 2 m s~* from use of the CORALIE spectrograph and 1.2-m Swiss telescope (ESO). 
De Ridder et al.l (120061 ) report ~3 m s _1 oscillations in the G9.5 III star e Oph. The power 
spectrum for the £ Hya radial velocities has maximum peaks corresponding to periods near 
0.15 days, and a generally complex structure with evidence for even mod e spacing expected 
of solar-like, coherent oscillations. The predicted mo de lifetimes from iHoudek fc Gough 



( 120021 ) are 17 days for thi s star, while IStello et al.l (120061 ) find a value of about two days from 
analysis of the extensive iFrandsen et al.l (120021 ) observations. The latter value, if true and 
common, could well limit the overall utility of asteroseismology of red giant stars. 

Photometric observations have been quite e f fectiv e for elucidating variations in extreme 



red giants to supergiants, e.g. iKiss Bedding! (120061 ) find ubiquitous variations in super- 



giants from multi-decade visual observations compiled by the AAVSO that are suggestive 
of stochastically dri ven oscillations. The O GLE-II data base with extensive data over three 
years also provides (IKiss fc Bedding! 120031 ) clear detection of multiple variability sequences, 
presumably related to different radial orders of oscillation for thousands of stars below the 
tip of the red giant branch from the LMC For these upper red giant branch variables ampli- 
tudes are commonly 1-4% with periods of 15-20 days with excitation mechanisms related to 
either Mi ra-like pulsators or sto c hastic ally-driven pulsations from convection. Farther down 
the RGB Edmonds fc Gillilandl (Il996h claimed that red giants often showed variations on 
time scales of ~2-4 days with amplitudes o f 0.5-1.5% for K gian ts in 47 Tuc based on a 40- 
hour U-band sequence obtained with HST. iJorissen et al.l (119971 ) used extensive Stromgren 
photometry of unusually high photometric stability to detect variations in K giants down to 
~0.5%, and, significantly, to establish that a clear level of var iability onset exists as a function 
of stellar parameters above which all red giants are variable. iHenry et al.l (120001 ) reported on 
extensive auto matic telescope observ ations of 187 red giants with sensitivity to about 0.2% 
confirming the IJorissen et al.l (119971 ) results and further establishing that most red giants 
earlier than G2 an d later than K2 are variables. At the 2 mmag precision level available 
Henry et al.l (120001 ) found that most red giants over G3-K1 were not variable. The variables 
to the red were interpreted as radial pulsations while those in the blue set had characteristic 
time scales too long for radial pulsations and non-radial g-modes were proposed. 

Extensive ground-based campaigns attempting to detect oscillations in red giants over 
luminosity, r adius ranges covered by the current data have largely fallen short of success. 
Stello et al.l (120071 ) observed several gi ants in M67 and report ed indications of excess power 
broadly consistent with expectations. IFrandsen et al.l (120071 ) observed many giants in the 
globular cluster M4 and reported no consistent evidence of oscillations. 



Given the extreme range of stellar properties pertaining to red giants, it is not surpris- 
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ing that multiple physical mechanisms may be responsible for subsets of these. With deep 
surface convection zones red giants are expected to show stochastically excited solar-like pu l- 
sations (IChristensen-Dalsgaard fe Frandsenl Il983l ; iHoudek et al.l Il999l ; ISamadi et al.l 120071 ) . 
Xiong fe Dens (120071 ) argue that the hotter type G giants fall within the classical Cepheid- 
S Scuti instability strip with k mechanism driving, and that a second instability strip for the 
cooler K and M giants at logfTsff) < 3.7 is driven by coupling between c onvection an d oscil- 
lations. The Kiong fe Deng (120071 ) result is broadly consistent with the iHenry et al.l (120001 ) 
observational results in finding that an intermediate temperature domain is pulsationally 
stable. 



The promise of space-based photometry has recently been demonstrated by lStello et al. 



(120081 ) with the detection of power spectrum excess in 11 red giants observed with the WIRE 
satellite. The WIRE observations of 15-61 days reach noise levels in amplitude spectra to 
a best of 7 ppm with a median of just over 20 ppm. Even more impressive are results 
from the MOST satellite on e Oph for which a 28 day observation set for this red giant 
with L/L Q ~ 60, T c ff ~ 4900 provides clear and unambiguous detectio n of multiple, equally 
spac ed modes from which s trong asteroseismic constraints are quoted by lBarban et al.l (120071 ) 
and iKallinger et al.l (120081 ) . Arguments for relatively long mode lifetimes of 10-20 days are 
encouraging. 

The data available for the present analyses were obtained with HST ACS as a nearly 
continuous time series of alternating F606W and F814W exposures for seven days in February 
2004. These data were collected for the purpose of e xtrasolar planet detection for Hot 
Jupiters orbiting upper main sequence stars in the bulge (jSahu et al.ll2006l ). The results here 
will be of particular interest for extending detections to smaller amplitudes than previously 
possible from ground-based observations, and for a much larger sample than possible with 
previous space-based results. 

Section 2 will provide an introduction to the properties of these HST data and discuss the 
techniques used for time series extractions on saturated stars. The general noise properties of 
stars on the RGB and main-sequence will be compared at comparable brightness levels and 
used to show in §3 that variations in the red giants are ubiquitous at some magnitude levels. 
The distribution of oscillation characteristics over the CMD will be developed in §4. Section 5 
will provide mappings of RGB position on the CMD to physical stellar parameters of mass, 
radius, luminosity and temperature and comparison of theoretically expected oscillation 
properties with those observed. Section 5 will also quantify how the variations detected 
along the galactic bulge RGB scale with L/M, compare with theory, and use this to predict 
expected variations that will be detectable by the Kepler Mission. 
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2. HST Time-Series of Galactic Bulge Stars 

The GO-9750 program from which these time series have been obtained was executed 
over 2004 February 23-29 using the ACS/WFC in filters F606W and F814W. The exposure 
times were always 339 seconds (the minimum providing high duty cycles with the WFC) and 
supported mmag precision near the bulge tu rnoff as needed f or detection of 1-2% transit 



signals which were the primary application (ISahu et al.l 120061 ) . Use of the two filters was 



usually equally split during the 105 HST orbits. The guiding accuracy of HST is excellent, 
but it can experience drifts on a timescale of a day comparable to the 0'/05 scale of the 
WFC pixels. To avoid potential systematic drifts in the photometry following from slowly 
developing pointing errors the observations were dithered each orbit by a random value 
intended to scramble offsets in time which would then allow derivation of intensity versus 
offset dependence and subsequent correction for this. The dithers plus inherent drift resulted 
in offsets that densely sampled the scale of a single pixel without exceeding this — exactly 
as desired. A total of 254 exposures were acquired in F606W and 265 in F814W providing 
an average of about five samples per 96-minute orbit during the 50 minutes available for 
science. 



Full difference image analysis (e.g. as used previously by iGilliland et al.l 120001 for a 



comparable WFPC2 project) was utilized for the bulk of the ti me series work. This involves 



build ing up a model for the intensity each pixel would record ( iGilliland. Nugent fc Phillips 



19991 ) under the assumption it is monitoring a non- variable target as a function of x, y offsets 
image-to-image and focus changes. Derivation of the image-to-image offsets and relative 
focus variations, as well as elimination of cosmic rays is an inherent part of the difference 
image analysis (DIA). However, DIA only works for unsaturated objects. Large step function 
gradients in flux associated with bleeding from saturated targets are not amenable to such 
modeling, and I reverted to a simple aperture sum in the direct images for the saturated stars. 
Apertures were generated that encompassed all of the pixels bled into from the saturation, 
the most critical factor in maintaining good results, but with an odd shape tracing the 
maximum extent of bleeding for a given star with a buffer all around of at least one pixel. 
Typical aperture sizes ranged from about 50 pixels for mildly saturated stars to over 600 
pixels at 5 magnitudes beyond saturation onset. 

All analyses were initiated with the pipeline provided _flt.fits files that contained data 
already corrected for contemporaneous biases and darks, and that had been flat fielded and 
converte d to electrons d etected. Using repeated exposures obtained for an HST calibration 



program 



Gillilandl (120041 ) showed that the photometric response for the ACS WFC remained 
perfectly linear (to precisions near the 0.1% level allowed by the data) up to and well beyond 
saturation when the GAIN = 2 option had been selected, which allows the device's full well 
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depth to be reached as the limiting factor rather than the analog to digital converter. One 
consequence of using simple aperture sums for the saturated stars rather than full DIA is 
that cosmic rays cannot be eliminated at the individual pixel level. Since a typical cosmic 
ray has 2000 electrons, cosmic rays can be eliminated as several sigma deviations even on 
pixels near the full well depth of 84,000 electrons with a good DIA model. However, typical 
saturated stars may have 4 x 10 6 electrons for which a 2,000 electron cosmic ray is exactly a 
1 a perturbation — not enough to be flagged and eliminated, but given that the typical star 
will experience an average of about one per exposure, enough to significantly raise the noise 
level. Unusually strong cosmic rays can of course still be flagged and eliminated even in the 
integrated sums. 

Input to the time series production includes th e pipeline da ta and a master star list 



(ISahu et al.ll2006l ) formed from separate, DAOPHOT (IStetsonlll992l ) analyses on over-sampled, 
combined images in both F606W and F814W. For the red giants, shorter 20 second expo- 
sures obtained as part of GO-9750 were also analyzed to provide master list photometry and 
star positions for stars that would be saturated and thus sometimes blended in the standard 
339-second images. Photometry is provided in the Vega magnitude HST system for the ACS 
F606W (Broa d V) a nd F814W (I) system using the prescription and zeropoints as given by 



Sirianni et al.l (120051 ). Figure 1 shows the CMD for 230,000 stars derived from these deep, 
combined ACS images. Figure 2 shows a blow-up of the brightest few magnitudes of the 
CMD with superposed fiducial boxes that select subsets of stars by magnitude and color. 

Relative time series are formed using the summed counts in individual images minus 
the mean over all images and then normalized by this mean. Time series values are left in 
these direct relative units rather than being transformed to magnitudes. A simple measure 
of time series quality for the full ensemble of stars analyzed is provided by plotting the 
rms, or standard deviation against magnitude as shown in Figure 3. For the stars that do 
not saturate, difference image analysis provides results that closely approximate limiting 
precisions expected from a combination of Poisson fluctuations on the object and sky, plus 
readout noise over nominal apertures. As discussed earlier the saturated stars, those brighter 
than F606W ~ 18.72, or F814W ~ 17.73, have an intrinsic noise level higher by about 20%, 
consistent with the inability to eliminate most cosmic rays from their photometric sums. 
After setting the sums over specially tailored apertures that encompass all of the pixels bled 
into for saturated targets, plus a small buffer, the relative time series are then decorrelated 
(multiple linear regression) with vectors of x, and y offsets over time, a measure of focus 
changes, the sky background level, a goodness of fit parameter from the difference image 
analysis steps, and an ensemble sum over 50 bright, non- variable stars. 

The distribution of rms values for the brightest stars in Figure 3 show positive offsets 
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from limiting values from two effects: the already mentioned inability to remove cosmic rays 
that leads to a jump near F606W = 18.7, and intrinsic variability of red giants, especially 
for those at the brightest magnitudes. Justifying the last point will be the primary topic of 
the next section. 



3. Ubiquity of Variations in the Red Giants 

Fortunately for the purposes of establishing the reality of red giant variations in these 
data, given that the saturated objects clearly have excess noise as evidenced by the jump in 
rms distribution in Figure 3, the CMD contains a large number of foreground dwarf stars at 
the same intensity levels. 

Figures 4 through 8 will contrast photometry for one foreground dwarf and one red 
giant (each chosen at random) from the intensity level boxes corresponding respectively to 
F606W = 15.5, 16.0, 16.5, 17.0 and 17.5 as shown in Figure 2. Each figure has the same 
form with direct ti me series plott ed in panels for the dwarf and giant separately at the top, 



and Lomb-Scargle (IScargld Il982f ) periodograms (square root of this calibrated in units of 
parts per thousand, or if multiplied by 1.085, are mmag for this amplitude spectrum) at the 
bottom. The time series are continuous over a seven day period except for losing half of 
each consecutive 96 minutes to Earth occultation, and cover HJD 2453058.52 to 2453065.47. 
The time series of the foreground dwarfs are usually closely approximated as white noise. 
By contrast the time series of red giants at the same brightness level show greatly enhanced 
variations at low temporal frequencies. The amplitude spectra are plotted from 1.65 to 100 
/iHz, or correspondingly periods of seven days to 2.7 hours. Both the dwarf and red giant 
time series have no additional structure (not shown) out to the Nyquist limit at about 17 
minutes except for expected aliases near the HST orbital frequency of 174 /xHz. Variations in 
the red giants are characterized by slow drifts with similar amplitudes in F606W and F814W 
and cannot usually be represented as a single sinusoid. At 18th magnitude in F606W the 
red giant oscillations are no longer noticeable from such comparisons in these data. 

To assess the generality of detectable oscillations among the red giants studied here 
it is informative to examine a histogram of rms levels for one of the magnitude bins of 
Figure 2. Figure 9 shows the rms distribution at F606W = 16.5 for dwarfs and red giants 
summed over bins B and C of Figure 2 centered on the RGB. The distribution for the 
dwarfs is strongly peaked for a bin from 0.0008 to 0.0009 with only a modest asymmetry to 
the high side. By contrast there are essentially no red giants at this magnitude and color 
(3 of 161) with an rms <0.0009, and the peak in the distribution function for red giants 
is at an rms within 0.0011 to 0.0013. A lower limit to time series variability (in units of 
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rms) at any magnitude may be formed by taking the square root of the difference between 
the mean red giant and mean dwarf variance (square of rms). For example, averaged over 
all the stars shown in Figure 9 the red giants have an excess rms relative to the dwarfs of 
0.00074 (ignoring one giant and two dwarfs that have large rms values more than 0.0015 
above the median value). Means for the dwarfs and giants in Figure 9 are 0.00092 and 
0.00118 respectively. If a significant component of the dwarf variability is intrinsic, rather 
than just noise from the observations and analysis, then the variability ascribed to the red 
giants should be increased accordingly. At the levels of nearly one part in a thousand being 
discussed here the values for the dwarfs are likely to be dominated by observational noise 
rather than intrinsic oscillation. Therefore the lower limit to red giant variations established 
in the simple quadrature difference approach is probably a close approximation to actual 
levels. 



Jorissen et al.l (119971 ) and iHenry et al.l (120001 ) had argued that most red giants in re 



stricted regions of the CMD down to sensitivity levels of 0.5% and 0.2% respectively are 
variable. This survey has shown that red giants significantly farther down the red giant 
branch are typically variable to sensitivity levels reaching a factor of 10 smaller. 



4. Distribution of Variation Properties Over the CMD 



Simple physical arguments, e.g. as developed by iKjeldsen fc Bedding! (119951 ). lead to 
predictions (albeit poorly constrained at this point) of how oscillations should vary with 
stellar properties. In this section I develop the empirical variations for red giant oscillations 
with position in the CMD for F606W amplitudes, the frequency of peak power, the relative 
amplitudes of F606W and F814W variations, a measure of mode multiplicity, and a mea- 
sure of the large splitting (expectation of equal separations in frequency between successive 
p- modes). In the next section I will establish the stellar properties as the same function 
of position in the CMD used for collation of oscillation properties and briefly compare ex- 
pectations from theory to the observations of oscillation properties as a function of stellar 
parameters. 



4.1. F606W amplitudes 

In §3 I gave a definition for lower limits to red giant oscillation amplitudes measured 
in units of rms by forming the quadrature difference of means for giants and dwarfs at the 
same magnitude and assuming the dwarfs are non-variable. The results of extending this to 
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six magnitude by five temperature bins defined in Figure 2 spanning the red giant branch is 
shown in Table 1. The excess rms is shown in units of parts-per-thousand. The equivalent 
amplitude of a single sinusoid yielding this variation level would be \/2 times the rms excess. 
Quantities in parentheses in Table 1 provide a measure of the error on the mean for each bin. 
The table entries are the means of two independent comparisons of rms values, one based on 
means evaluated using values clipped to be within the median plus 0.0015, and one simply 
based on median values for both the main sequence and red giant stars. Usually the values 
based on medians and means are very close in value. 

A number of things are immediately obvious from Table 1: (1) At and above F606W = 
17.5 through 15.5, variability is securely detected in all bins based on differences of variances, 
even via this rather crude approach. (2) Within each column (with two exceptions out of 21 
with the miss well within l-a) roughly corresponding to different temperatures, amplitudes 
increase monotonically to higher luminosity (assuming stars are at roughly equal distances 
as seems secure from the coherence of features such as the red giant clump seen in the CMD). 
(3) At the same magnitude or luminosity, with only one exception out of 20 steps (and there 
it is within one-half a on mutual error bars), amplitudes increase to lower temperature. 

The bins in Table 1 (and throughout this paper) without numerical entries were avoided 
for two different reasons. At the two faintest bins for column A there are a combination 
of not many stars and developing overlap with the foreground dwarf distribution. At the 
two faintest bins in column E there would be only one red giant per bin. The total number 
of dwarfs in the six bins over luminosity entering the statistics for Table 1 is 340, the total 
number of red giants distributed over 26 bins is 1034. 

4.2. Mean frequency of variations 

The distribution of peak power over the CMD is shown in Figure 10 where medians 
of amplitude spectra in ppt are given for the bins spanning both main sequence (taken as 
controls) and red giants. For Figure 10 the combined time series in F606W and F814W were 
analyzed simultaneously. The variation of amplitudes introduced in the previous section and 
shown in Table 1 are easily visualized in Figure 10. Noteworthy are the nearly flat median 
amplitude spectra for the main sequence stars with just a hint of turn up (possible l//noise) 
below frequencies of 5 /iHz, or periods greater than two days. 

Inspection of Figure 10 strongly suggests that with increasing luminosity or decreasing 
temperature the location of peak power moves to lower frequency. I have quantified this by 
defining a range of frequency in each star spanning the first and last locations with amplitude 
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greater than or equal to half the peak amplitude up to 70 fiHz and then computing the mean 
frequency over the domain weighted by the amplitude spectrum. Only stars for which the 
power peak a t frequency le ss than 70 fiHz exceed a significance limit of 0.001 for chance of 
false alarms (jScargle 1982 ) are included in the statistics. Within each row and column of 
Table 2 the changes of typical frequency at which oscillations peak occurs in the expected 
way. 



4.3. Relative amplitudes of F606W and F814W variations 

The variation of relative amplitudes with bandpass over the CMD is shown in Table 3. 
These are computed by Gaussian smoothing and interpolating each of the F606W and F814W 
time series to a common grid (with full width half maximum of Gaussian at 2 hours), then 
performing a linear fit of F606W to F814W. Stars are retained using three cuts: as for 
Table 2 and mean frequency the star must have a power spectrum peak with a false alarm 
probability less than 0.001 and the peak must be at less than 70 /iHz, in addition, the linear 
correlation of F606W and F814W filtered time series must be >0.6. Here distinct trends are 
clear with larger amplitudes in F606W relative to F814W at brighter and cooler stars. 



4.4. Mode multiplicity measure 

Most red giants that show any power spectrum peak with a false alarm probability of less 
than 0.001 show multiple peaks in the power spectrum. To quantify this I have counted the 
number of peaks below 70 /zHz at an amplitude greater than four times the mean amplitude 
spectrum averaged over 70-100 /d3z. This algorithm returned 3, 3, 4, 2, 1 peaks for the 
red giants plotted in Figures 4 through 8 respectively. None of the dwarf stars plotted in 
these figures showed any peaks above this mean level. Table 4 shows the distribution of 
mean number of peaks over the CMD. Generally lower values in the left-most column and 
faintest row follow from having barely enough S/N to detect oscillations. The number of 
independent peaks for stars toward the upper right corner of the table are limited by the 
frequency resolution of these data being insufficient to resolve independent peaks (see the 
pile up at lowest frequency for these bins in Figure 10). The only conclusion is that red 
giants that show variations are characterized by changes on several independent frequencies 
rather than commonly displaying only one or two modes. 
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4.5. Large splitting of modes 

A power spectrum of solar oscillations shows a well populated set of modes at nearly 
equal frequency separation. The most evident separation is called the "large splitting" and 
corresponds to modes that differ by one in radial number, or two in the number of nodes 
around the stellar circumference. Measuring the large splitting for stars displaying solar-like 
oscillations immediately gives a good constraint on the mean density of the star. 

Inspection of Figure 4 is suggestive of three modes at equal separation. I have attempted 
to quantify the large separation for the 200 red giants that each showed at least four peaks 
at greater than four times the mean amplitude spectrum averaged over 70-100 /iHz by 
computing the power spectrum of the power spectrum over the frequency domain spanning 
these significant peaks. Unfortunately for every red giant like the one in Figure 4 with 
regularity in the power spectrum, there are several more red giants for which the power 
spectrum of the power spectrum fails to reflect such regularity. A formal evaluation of 
median large splittings as the location of the highest peak in the power spectrum of the 
power spectrum up to 15 /iHz fails to show any trends over the CMD and I have elected not 
to show the resulting table. 

These time series are relatively short for quantifying frequency characteristics of the 
oscillations in red giants. Reliable detections of large splittings for a large ensemble of red 
giants (see discussion of prospects with Kepler in §5.3) awaits collection of more extensive 
observations. It remains to be seen whether red giants will typically have mode lifetimes 
sufficiently long to allow frequent measurement of the large splitting. 



5. Relation to Stellar Physical Properties 



In the previous section observational results were developed that show clear correlations 
of variation amplitudes, relative amplitudes in F606W and F814W, and mean frequency of 
variations for a wide range of galactic bulge red giants over the CMD. In order to provide 
interpretation of these measured variation/oscillation properties, and to make predictions for 
other characteristics for which these data were insufficient, it is necessary to have estimates 
of stellar physical parameters for the observed stars. Lacking any detailed knowledge for 
individual stars in this field beyond photometry in two HST ACS bandpasses, I make use of 
the easy to use and conceptually power ful simulation package B a STI (for A Bag of Stella r 
Tra cks and Isochrones) documented in iPietrinferni et al.l (120041 ). iPietrinferni et al.l (120061 ). 
and lCordier et all toQ^ \. 



In particular I have utilized the BaSTI web tool for generating synthetic color-magnitude 
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diagrams with input appropriate for the galactic bulge population. The option for simulating 
the M ilky Way Bulge star formation history was selected; this is based on lMolla. Ferrini &: Gozzi 
( 120001 ) . The synthetic CMD provided by BaSTI included only standard Johnson- Cousins 
photometric bandpass simulation. In order to transform to the ACS F606W and F814W 
system in use for these observations I ran identical isochrone simulations with BaSTI ap- 
propriate to the galactic bulge RGB for which output in both Johnson-Cousins and ACS 
bandpasses were available, and derived a simple empirical transform as n eeded to take the 
CMD simulations to the ACS system. Having input values of 14.3 and 0.64 (ISahu et al.ll2006l ) 
for distance modulus and E(B-V) color excess for the BaSTI CMD simulation, transforming 
photometric results to the ACS system a further tweak of 0.05 in F606W-F814W provided 
the simulated photometry shown in Figure 11 which compares well with the observations in 
Figure 2. Random errors of 0.1 magnitude in F606W and 0.03 in F606W-F814W were also 
applied to the BaSTI simulations shown in Figure 11. 



5.1. Mass, radius, luminosity, temperature over the CMD 

For the purposes here of developing a rough understanding of how intrinsic stellar prop- 
erties vary over this observational data set the matching simulation from BaSTI shown in 
Figure 11 allows easy derivation of mean physical parameters corresponding to the boxes in 
F606W versus F606W-F814W used throughout this paper. 

The output from BaSTI includes age, [Fe/H], luminosity, mass, effective temperature 
and apparent photometry for a large ensemble of simulated stars. Having transformed the 
photometry to the ACS system it is now a simple matter to compile averages within the 
CMD boxes for stellar physical properties of interest. 

Table 5 shows mean luminosities in solar units from the BaSTI simulations. Over the 
range of F606W = 17.5 in box B of the CMD through F606W = 15.5 in box E for which 
RGB variations are securely indicated for the majority of stars the stellar luminosities range 
from 15 to 350 solar, or a dynamic range of more than 20. 

Effective temperature means are shown in the middle block of Table 5 with a rang e from 



abou t 5100 K to 3950 K corresponding to spectral types of about G4 to K3 giants ((Allen 



19731). 



Stellar radii are shown in the lower block of Table 5 with a range from about 5 to 
40 in units of the solar radius. Not shown, but available from BaSTI are results for mass 
and intrinsic metallicity. The mass values vary little, being about 0.9 M & at the highest 
temperatures to about 1.05 M at the lowest temperatures. Metallicities range from —1.0 
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for box A at the left side of the RGB to 0.0 for box E at the right side. The [Fe/H] values are 
unique in showing relatively large distributions within these CMD boxes, with dispersions 
within boxes often near a factor of two. Other quantities like luminosities, temperatures, 
radii and masses are much more tightly coupled to CMD position with dispersions within 
boxes remaining small compared to mean changes box-to-box within the CMD. 



5.2. Expected oscillation properties compared to observations 



With values of typical stellar properties of luminosity, radius and effective temperature 
shown in the previous section the corresponding theoretically expected oscillation properties 
will be developed here. 



For the Sun, p-modes show a broad distribution centered on 3020 fiHz. iBrown et al. 



( 119911 ) argued that the location of maximum power would likely remain constant in ra- 
tio to the acoustic cut-off frequency for other stars. This may be expressed as (following 
Kieldsen fc Beddind[l995} ): 



3021- 



M/M G 



(R/RqY (T eff /5777) 



1/2 



jjiHz 



Using inputs from Table 5 plus mass values differing only slightly from one solar mass 
(but included) the predicted location of maximum power from Eq. (1) is shown in the top 
block of Table 6. Note that the variation from lower left to upper right in the observational 
version of this (Table 2) is in good agreement as to simple trends. In detail the theoretical 
values of Table 6 decline to much shorter frequencies for the most luminous and cool stars 
than is seen in the short observation sets available here. 

As expected theoretically, but still wi th considerable uncertain t y, the most likely mo des 
of oscillation may be low-order p-modes (jDziembowski et al.ll200ll ; Kiong fc Deng! 120071 ) in 
which case the f reque ncies would be somewhat different than given by Eq. (1) and Table 6. 
Kallinger et al.l (120051 ) presented results for modelling a 47 Tuc red giant with L/L Q = 119, 
R/R Q = 18.8, and M/M Q = 0.9 which very closely matches bo x C at F606W = 16.0. With a 



<4% adjustment for square root of relative density between the iKallinger et al.l (120051 ) model 
and this box I estimate the 3rd radial overtone frequency at 8.42 /iHz. Scaling this result by 
< p> x / 2 leads to the values shown in the middle block of Table 6. These frequencies are a 
considerably better match to the location of power maxima shown in Table 2. Allowing the 
freedom to select different radial overtones from the fundamental through 5th as a function 
of CMD position would allow essentially perfect matches. 
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The expected large splitting for modes that differ by one in radial order, or two in 
azimuthal is inversely proportion to the stellar mean density (with correction factors of only 
a few percent needed). The modes of even and odd azimuthal order produce an apparent 
evenly spaced set of modes separated as: 

Az/ /2 = (M/M ) 1/2 (R/R & )~ 3/2 67.45 /iHz (2) 



The lowest block in Table 6 shows the implied separation of equally spaced non-radial 
modes. The extent of the observations provides a resolution of 1.65 yuHz, the upper-right 
half of the table has large splittings that are smaller than the resolution available. While 
these observations are useful for broadly illuminating properties such as typical amplitudes, 
they are clearly not sufficient to pursue this interesting quantity required to perform basic 
asteroseismology. At the magnitude-temperature combinations in Table 6 where splittings 
exceed 3 /zHz and might be expected to be robustly detectable in these data from a resolution 
perspective, Table 1 shows that there is very little overlap with respect to clear indications 
of excess power. 

For low-order radial modes typical frequency separations would be 1.5-2 times those 
shown in Table 6. Over much of the CMD domain the resolution provided by these data 
would still be insufficient to separate individual modes. However, stars at F606W = 17.0- 
17.5 for boxes B and C might be amenable to more detailed study — a topic for future 
investigation for individually favorable stars. 

The remaining lowest order characteristic of oscillations worth comparing theoretical and 
realized results for is the variation amplitude. For this the literature indicates significantly 
less confidence in provid ing a unique prediction. I have adopted two predict ions in units of 
parts-per-thousand from lKjeldsen fc Bedding! (119951 ) and lSamadi et al.l (120071 ) with equations 
of: 

-3 L/Lq 



(SL/L) = 4.27 x 10" 



(T cff /5777f M/M(? 



(3) 



(SL/L) = 4.27 x 10~ 3 (L/M) 



0.7 



(4) 



The results from Eq. (3) of lKieldsen fc Bedding] (119951) are shown in the upper block of 
Table 7, while those from Eq. (4) from ISamadi et all ( 2007 ) are given in the lower Table 7 
block. Comparing these two tables to Table 1, which shows the excess rms of the observed 
time series, demonstrates excellent agreement in terms of trends over the CMD. Taking into 
account that rms excess is not the same to within factors of a few as the theoretically given 
oscillation amplitudes, it is not possible to strongly discriminate between the two theoretical 
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predictions. Eq. (3) generally suggests larger amplitudes than exist in the data. Eq. (4) 
lacks the strong dependence on T e g that seems to be present in the data, and has amplitudes 
that are too small. Furthermore both theoretical predictions show a relatively larger spread 
within columns over the magnitude range than is seen in the data, e.g. for column D over 17.5 
to 15.5 the observations imply a factor of 4.0 change while the two t heoretical evaluations 



yield 12.8 and 5.0. In the latter regard the simpler L/M scaling of the ISamadi et al.l (120071 ) 
relation is much closer to these observations, and the existence of many independent modes 
summing incoherently can explain at least part of the offset between absolute values. 



The general rise of amplitudes to 
consistent with observational claims of 



ower effective tem perature in these data is broadly 



giants, and with theory as in Kiong fc Deng! (120071 ). The distribution of red giants in this 



Henry et al.l (120001 ) for generally more luminous red 



study does not extend to early G where amplitudes might be expected to rise to higher 
effective temperatures. 

The observations here for which excess power has been detected in most time series 
span a range of 20 in L/M, and it will therefore be informative to attempt to develop 
fits for the amplitude as a function of stellar parameters from these data. I ran experiments 
performing linear regressions of the rms excesses of Table 1 against logarithmic combinations 
of luminosity, mass, temperature and radius in various combinations. The reduced x 2 versus 
a fit in just L/M is 7.74 — consistent with the large scatter in the left panel of Figure 12. 
From inspection of Figure 12 the amplitude excesses show systematic offsets corresponding 
to left-right position within the red giant branch. Since left-right position is associated with 
strong temperature changes a regression in L/M and T e g/ '5777 was performed dropping the 
reduced x 2 to 1.68 as shown in the right panel of Figure 12. This formal fit is: 

Amplitude = -1.34 + 0.483 log 10 (L/M) - 14.2 log 10 (T ofr /5777) (5) 

Since the amplitude fit here is only a rough indication of excess variation from all sources, 
not just potential coherent oscillations as might be revealed by longer and more precise 
observations, general adoption of this equation is not intended to represent a prediction of 
coherent oscillation amplitudes in red giants. 



5.3. Significance for Kepler Mission 

These observations have been powerful in indicating the ubiquity of variations with 
characteristic time scales similar to those expected from stellar oscillations in nearly all red 
giant stars observed with luminosities above 20-30 L/L Q . In terms of providing results of 
direct asteroseismic significance, e.g. through determination of the large splitting for solar-like 
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p-mode oscillations which would yield a constraint on stellar density, these results are severely 
limited by both the precision and frequency resolution which are generally inadequate to 
provide such results. 



The Kepler Mission (IBorucki et al.l 120071 ) which has a primary goal of detecting Earth- 



like planets transiting roughly solar-like stars also carries a goal of characterizing the host 
stars in order to better determine planet parameters and the context in which they are 
formed and exist. Two key tools to be used by Kepler in characterizing a large subset of 
the 100,000 stars to be observed nearly continuously for 3.5 years will be the derivation of 
astrometric parallaxes, and asteroseismic probing of interiors. Ninety-nine percent of the 
stars to be observed will be chosen on the basis of yielding the best possible chances for 
detecting small orbiting planets, thus restricting such choices to dwarfs. However, current 
plans call for observing 1,000 red giants at the standard 30-minute cadence for 3.5 years 
for the dual purposes of: (a) providing low parallax controls, and (b) asteroseismic results 
expected to much better constrain our knowledge of these stars and thus contribute to a 
better overall understanding of stellar characteristics within the Kepler field. In order to 
serve as astrometric standards it is important that targets not be allowed to saturate in the 
individual 3-8 second (~6 seconds currently favored) exposures which will be co-added to 
30 minutes on board. This sets a magnitude limit of about 11.7. A test selection of 1,000 
red giants distributed over the focal plane chosen to never saturate, weighted to be at large 
distance, and generally lacking in close neighbors that would compromise astrometry have a 
mean magnitude of 11.98 (in a system appropriate to the broad bandpass of Kepler from 400 
to 850 nm). At this magnitude the expected precision will be 65 ppm per 30 minute sum. 
After a one- week period the Kepler photometry should yield power spectrum noise levels of 
3.5 ppm for these 1,000 red giants compared to a best noise level of about 35 ppm in these 
HST observations of one week. More importantly the Kepler observations are expected to 
be nearly continuous for a full 3.5 years. After the acquisition of three months of data the 
precisions provided for stars near 12th magnitude should reach about 1.0 ppm (sufficient 
to detect oscillations in solar-like stars, let alone the red giant variations at much larger 
amplitudes), and provide a frequency resolution adequate to critically sample the evenly 
spaced p-mode distribution for stars corresponding to all of the CMD positions detailed 
in this paper. Out of the 1000 red giants selected on criteria primarily intended to allow 
excellent astrometric measures from the Kepler data, application of Eq. (5) yields only 31 
stars that have predicted rms excesses below 100 ppm; the majority are expected to have 
excesses greater than 1000 ppm. 

With a launch expected in February 2009 Kepler should quickly provide a vast wealth 
of information bringing clarity to properties of oscillations in red giants that have only been 
weakly explored to date. 
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6. Summary 

Observations over a seven day span with the ACS imager on HST have been used to 
demonstrate that most red giants with luminosities between 30 and 350 L Q show evidence 
of excess noise in power spectra in comparison to main sequence control stars at the same 
brightness. Amplitudes and frequencies depend upon position in the CMD, and thus differ- 
ences of stellar parameters of L, M and T c g- in rough agreement with theoretical expectations. 
The best match of frequency for observed power excess in the red giants is consistent with 
interpretation as low-order radial modes. Since the frequency resolution of these data are not 
sufficient to resolve expected oscillation characteristics over much of the sampled domain, 
asteroseismic applications are not provided. The upcoming Kepler Mission will provide im- 
provements of more than an order of magnitude in both photometric precisions and frequency 
resolution allowing definitive understanding of the nature of variations in red giants all along 
the branch. 
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Fig. 1. — CMD for 230,000 stars identified from combined, over-sampled images consisting 
of 86,106 (F606W) and 89,835 (F814W) seconds with ACS on the galactic bulge. 
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Fig. 2. — Region of CMD from Figure 1 replotted on a larger scale to emphasize the red 
giant domain. The over-plotted grids will be referred to throughout the paper to demonstrate 
changing properties of variations as a function of stellar parameters. 
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Fig. 3. — Shows standard deviation of relative time series intensity for all 230,000 stars 
shown in Figure 1. The break near F606W = 18.5 results from transition from difference 
image analysis (fainter) to simple aperture photometry as discussed in the text. The solid 
curve shows expected limiting precision (Poisson on objects and sky plus readout noise) 
which closely tracks results for the non-saturated star photometry. 
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Fig. 4. — Compares time series for a foreground dwarf and red giant star at comparable 
magnitude of F606W ~15.5. Top two panels for the dwarf and red giant respectively show the 
direct, relative photometry spanning 7 days with results for bot h filters shown . The bottom 
two panels show square root of the Lomb-Scargle periodogram ( ]Scarglelll982l ) calibrated to 
units of parts-per-thousand. Note that the red giant shows excess variability over a domain 
of 2.5-10 /iHz (about 1-5 days). 
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Fig. 5.— Same as Figure 4, but for F606W ~16.0. 
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Fig. 6. — Same as Figure 4, but for F606W ~16.5. Excess variability extends dows to periods 
of about onedialf day. 
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Fig. 7. — Same as Figure 4, but for F606W ~17.0. Although excess variability is still 
nocticeable in the red giant, the amplitude is considerably smaller than is typical of giants 
at higher luminosity. 
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Fig. 8.— Same as Figure 4, but for F606W ~17.5. 
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16.5 and red giants at the same magnitude pulled from bins B and C. 
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Fig. 10. — Individual boxes map to the distribution shown in Figure 2 and Table 1, the 
frequency range and amplitude spectrum range are the same for each as shown for lower left 
case. The number of stars contributing to the median amplitude spectrum is shown in the 
upper right of each box. The left column is for the dwarf star control sample. 
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Fig. 11. — Same as Figure 2, but with simulated stellar properties from the BaSTI codes. 
See text for details. 
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Fig. 12. — Symbols and associated error bars correspond to the variation excess rms in parts 
per thousand as given in Table 1. In the left panel these are plotted against the luminosity 
to mass (in solar units) ratio developed from the analyses discussed in §5.1. The right panel 
plots the same data against a least squares regression in L/M and T e g as discussed in the 
text and developed in functional form as Eq. (5). 
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Table 1. Excess rms of Red Giants dependence on CMD position. 



F606W 


A 


B 


C 


D 


E 


15.5 


0.64(.16) 


0.91(.26) 


1.07(.18) 


1.92(.24) 


2.59(.37) 


16.0 


0.59(.09) 


0.87(.09) 


1.27(.12) 


1.79(.12) 


2.69(.25) 


16.5 


0.54(.13) 


0.68(.12) 


0.80(.12) 


1.39(.15) 


1.90(.ll) 


17.0 


0.39(.17) 


0.59(.10) 


0.76(.08) 


1.01(.18) 


1.34(.14) 


17.5 




0.41(.10) 


0.54(.08) 


0.48(.12) 




18.0 




-0.11(.20) 


-0.10(.15) 


0.71(.21) 





Note. - See Figure 2 for correspondence of bins in table with 
position in the CMD. From A to E spans lower to higher F606W- 
F814W color, or higher to lower temperatures. Units are parts-per- 
thousand of measured rms excess of giants relative to dwarfs at the 
same F606W, errors on the mean in each bin are shown in parentheses. 



Table 2. Median of mean frequency - dependence on CMD position. 



F606W 


A 


B 


C 


D 


E 


15.5 


17.7(5) 


11.05(4) 


7.2(8) 


4.0(6) 


5.0(2) 


16.0 


15.4(21) 


12.2(21) 


8.4(10) 


7.2(9) 


5.2(10) 


16.5 


18.9(43) 


17.5(70) 


16.3(43) 


9.2(18) 


6.7(5) 


17.0 


25.2(9) 


18.4(41) 


19.2(58) 


16.6(12) 




17.5 




31.6(4) 


23.3(17) 


25.3(2) 





Note. — Values of the median frequency in amplitude spec- 
tra for each CMD bin with at least 2 star s having signi ficance 



of highest peak up to 70 /iHz < 0.001 flScardd 119821 ) . The 
number of contributing stars is shown in parenthesis. Units 
are /xHz. 
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Table 3. Relative amplitude of F606W to F814W - dependence on CMD position. 



F606W 


A 


B 


C 


D 


E 


15.5 


0.96(0.12,3) 


0.99(0.19,4) 


1.14(0.23,8) 


1.40(0.42,6) 


1.65(0.18,2) 


16.0 


0.97(0.18,15) 


1.07(0.12,20) 


1.33(0.16,10) 


1.37(0.24,9) 


1.72(0.22,8) 


16.5 


0.93(0.19,30) 


1.01(0.18,61) 


1.11(0.25,37) 


1.44(0.20,18) 


1.73(0.23,5) 


17.0 


1.14(0.10,5) 


0.96(0.19,33) 


1.08(0.19,48) 


1.14(0.19,11) 




17.5 






0.96(0.16,10) 


1.07(0.05,2) 





Note. - Relative amplitude of F606W time series variations to same in F814W, 
see text for details. Standard deviation and number of stars included in each bin are 
shown in parentheses. 



Table 4. Mode multiplicity - dependence on CMD position. 



F606W 


A 


B 


C 


D 


E 


15.5 


4.40 


5.00 


3.75 


2.50 


4.50 


16.0 


3.48 


4.45 


5.00 


4.22 


3.40 


16.5 


2.77 


3.53 


4.35 


4.00 


4.60 


17.0 


2.78 


2.66 


3.12 


4.33 


5.00 


17.5 




1.25 


1.82 


1.00 





Note. - Mean number of significant 
peaks per power spectrum — see text for de- 
tails. 
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Table 5. Stellar Characteristics over CMD - BaSTI simulation. 



F606W A B C D 


E 


Mean Luminosity (solar units) 


15.5 


164.6 190.1 217.8 270.4 


350.3 


16.0 


84.2 109.1 123.3 148.6 


277.3 


16.5 


57.7 57.1 58.2 84.5 




17.0 


45.7 51.1 52.5 55.4 




17.5 


25.9 25.0 27.2 28.8 




18.0 


14.7 15.2 16.6 21.8 




Effective Temperature (Kelvins) 


15.5 


4887 4550 4385 4093 


3949 


16.0 


4943 4677 4457 4246 


3999 


16.5 


5176 4819 4688 4375 




17.0 


5000 4819 4699 4508 




17.5 


5058 4977 4710 4667 




18.0 


5135 5093 4819 4732 




Mean Stellar Radii (solar units) 


15.5 


17.97 22.28 25.66 32.83 


40.14 


16.0 


12.56 15.97 18.70 22.61 


34.82 


16.5 


9.48 10.88 11.61 16.06 




17.0 


9.04 10.29 10.97 12.25 




17.5 


6.65 6.75 7.87 8.24 




18.0 


4.86 5.02 5.86 6.98 





Note. — See text for discussion of simulations 
based on BaSTI codes. 
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Table 6. Predicted frequency characteristics. 



F606W 


A 


B 


C 


D 


E 


Frequency of maximum power (/-iHz) 


15.5 


7.83 


5.89 


4.84 


3.35 


2.36 


16.0 


17.48 


11.49 


9.39 


7.07 


3.13 


16.5 


28.00 


25.45 


24.34 


14.23 




17.0 


34.27 


28.88 


27.24 


23.81 




17.5 


63.05 


63.35 


55.71 


52.39 




18.0 


117.53 


112.16 


99.17 


71.89 




Frequency of 3rd radial overtone (/iHz) 


15.5 


8.11 


6.13 


5.20 


3.74 


2.80 


16.0 


14.45 


10.29 


8.42 


6.65 


3.54 


16.5 


21.32 


18.61 


17.57 


11.23 




17.0 


23.92 


20.38 


19.13 


16.74 




17.5 


38.06 


37.64 


32.23 


30.46 




18.0 


60.93 


58.44 


50.22 


38.99 




Large frequency splitting (/iHz) 


15.5 


.78 


.59 


.50 


.36 


.27 


16.0 


1.39 


.99 


.81 


.64 


.34 


16.5 


2.05 


1.79 


1.69 


1.08 




17.0 


2.30 


1.96 


1.84 


1.61 




17.5 


3.66 


3.62 


3.10 


2.93 




18.0 


5.86 


5.62 


4.83 


3.75 




Note. 


- All quantities 


use the 


BaSTI 


derived 



stellar parameters as shown in Table 5. Evalua- 
tion of the predicted location of maximum oscillation 
power uses Eq. (1). S caling of 3rd overtone frequency 
( Kallinger et al.fl"2005[ ) by p 1 / 2 provides an indication 
of low-order radial mode frequencies. Prediction for 
separation of modes in frequency space as half of the 
large splitting is from Eq. (2). 
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Table 7. Predicted oscillation amplitudes. 



F606W 


A 


B 


C 


D E 


Kieldsen & Bedding 


(1995) 


- Eq. (3) 


15.5 


1.28 


1.52 


1.76 


2.29 3.07 


16.0 


.58 


.81 


.93 


1.15 2.36 


16.5 


.39 


.38 


.39 


.59 ••• 


17.0 


.30 


.34 


.35 


.37 ••• 


17.5 


.17 


.16 


.17 


.18 •■■ 


18.0 


.09 


.09 


.10 


.13 •■■ 


Samadi et al. (2007) - Ea. (A) 


15.5 


.18 


.19 


.20 


.21 .25 


16.0 


.11 


.13 


.13 


.14 .21 


16.5 


.09 


.08 


.07 


.09 


17.0 


.07 


.07 


.07 


.07 •■■ 


17.5 


.05 


.04 


.04 


.04 


18.0 


.03 


.03 


.03 


.04 • • • 



Note. Amplitudes in parts-per- 

thousand for comparison with observed 
rms given in Table 1. 



